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ABSTRACT 



We discuss how different tfieoreticai predictions for the variance cr^ of the counts- 
in-cells distribution of radio sources can be matched to measurements from the FIRST 
survey at different flux hmits. The predictions are given by the integration of models 
for the angular correlation function w{6) for three different functional forms of the 
redshift distribution N{z), different spatial correlation functions that match the ob- 
served present day shape and by different evolutions of the bias b{z) with redshift. We 
also consider the two cases of open and flat Universes. Although the predicted w{9) 
show substantial differences due to differences in the A^(z)'s, these differences are not 
signiflcant compared to the uncertainties in the current observations. It turns out that, 
independent of the geometry of the universe and the flux limit, the best fit is provided 
by models with constant biasing at all times, although the difference between models 
with epoch-independent bias and models with bias that evolves linearly with redshift 
is not very large. All models with strong evolution of bias with epoch are ruled out 
as they grossly overestimate the amplitude of the variance over the whole range of 
angular scales sampled by the counts-in-cells analysis. As a further step we directly 
calculated Wohs(^) at 3mJy from the catalogue and matched it with our models for 
the angular correlation function in the hypothesis that the clustering signal comes 
from two different populations, namely AGN-powered sources and starbursting galax- 
ies. The results are consistent with a scenario for hierarchical clustering where the 
fainter starbursting galaxies trace the mass at all epochs, while the brighter AGN's 
are strongly biased, with h{z) evolving linearly with redshift, as suggested by some 
theories of galaxy formation and evolution. 

Key words: galaxies: clustering - radio galaxies - cosmology: theory - large-scale 
structure 



1 INTRODUCTION 

During the past few years there has been an increasing inter- 
est in how galaxies trace mass and specifically in the issue 
of a biasing factor evolving with redshift (e.g. Fry, 1996; 
Bagla, 1997; Matarrese et al., 1997; Moscardini et al., 1997; 
Tegmark & Peebles, 1998). This has received added impe- 
tus after Steidel et al. (1996, 1997) reported evidence for 
a strong concentration of galaxies at z ~ 3 (Lyman-Break 
Galaxies) , which would imply a value of 6 ~ 6 at those red- 
shifts. 

Radio objects can be detected up to very significant cosmo- 
logical distances {z ~ 4) and therefore provide information 
on large scale structure at rather early epochs when the main 
growth of structures occurred. Even though evidence of clus- 



tering in radio catalogues was detected in early wide-area 
surveys (Seldner & Peebles, 1981; Shaver & Pierre, 1989; 
Kooiman et al., 1995; Sicotte, 1995; Loan, Wall & Lahav, 
1997), the FIRST survey (Becker et al., 1995) is the first one 
in which the angular clustering of radio sources down to the 
mjy level is detected with high signal to noise ratio (Cress et 
al., 1997; Magliocchetti et al., 1998). Unfortunately the rela- 
tion between angular measurements and spatial quantities is 
strongly dependent on the radio source redshift distribution 
N{z), which becomes more and more uncertain as the flux 
threshold is lowered. Dunlop and Peacock (DP, 1990) pro- 
vided models of epoch-dependent luminosity functions for 
radio sources to make estimates of N{z) that work well for 
bright sources, but these predictions drastically diverge at 
faint fluxes. One of the main uncertainties at such low flux 
densities is given by the presence and relative amplitude of a 
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low-redshift spike due to a population of starbursting galax- 
ies (e.g. Windhorst et al., 1985) now believed to constitute 
a majority of sources at mjy levels (see Wall, 1994 for an 
overview) . 

The knowledge of a well assessed functional form for the red- 
shift distribution of radio sources N{z) has become of cru- 
cial importance in the last decade for both radio-astronomy 
and cosmology. In the former case this would provide tests 
for radio-source unification and luminosity evolution models 
(Wall & Jackson, 1997), while in the latter case it would al- 
low conversion of angular clustering measurements to spatial 
clustering estimates that can be used to constrain structure 
formation models. 

We present here a theoretical approach able to put con- 
straints on both the functional form of A^(2:) at fluxes 
F < lOmJy and the clustering properties of radio objects. 
We pay particular attention to the evolution of bias, i.e. the 
way radio galaxies trace the mass distribution. We begin by 
assuming that the spatial correlation function of mass is de- 
termined by the linear growth of mass fluctuations. This is 
a reasonable assumption because the mean redshift of radio 
sources in the survey (z ~ 1 as opposed to z ~ 0.1 obtained 
for optical and infrared surveys), and so the angular scales 
that we consider correspond to physical scales where the 
linear theory still holds. The clustering of radio galaxies is 
then related to mass via three different models for the evo- 
lution of bias with redshift, each one related to a reasonable 
model for galaxy evolution, as taken from Matarrese et al., 
1997 and Moscardini et al., 1997. We then predict the an- 
gular two-point correlation function w{9) for three models 
of N{z) that span the range realistic distributions for faint 
radio objects. We obtain the correlation function using both 
flat and open geometry. 

In section 2 we introduce the different models for N{z) used 
in the projection analysis, while section 3 describes the cal- 
culation of the theoretical predictions for w{9); the predic- 
tions for each model are then presented in section 4. In sec- 
tion 5 we move to the analysis of the variance and we show 
the results obtained for the measurements of the of the 
distribution of sources in the FIRST survey at flux limits up 
to 10 mJy. The data are compared with our predictions from 
different models of w{9) in section 6. Section 7 presents our 
results on the interpretation of clustering and bias under the 
assumption that radio sources are formed by two different 
populations. Section 8 summarises our conclusions. 



2 THE REDSHIFT DISTRIBUTION N(Z) 

To test the stability of the predicted angular correlation 
function given for different functional forms of the redshift 
distribution N{z), we use the Dunlop & Peacock (DP, 1990) 
models of epoch-dependent space density for radio sources. 
These models used Maximum Entropy analysis to determine 
polynomial approximations to the luminosity function and 
its evolution with redshift; the approach incorporated the 
then- available identification and redshift data for complete 
samples from radio surveys at several frequencies. 
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Figure 1. Redshift distribution N{z) for the radio source popu- 
lation at 1.4 GHz at a flux limit of 3 mJy. The dotted and dashed 
curves represent the 6 models (1-4, 6-7) taken from Dunlop & 
Peacock (1990); the solid curve is the average. Model 5 is omitted 
as it shows a totally unrealistic sharp and dominant feature at 

2 = 4 

Figures ^ and |^ show, for flux limits of 3 and 10 mJy 
respectively, models 1-4 and 6-7 taken from DP (dotted- 
dashed lines) and their average (solid line). The discrepan- 
cies amongst different models become increasingly signifi- 
cant as the ffux threshold is lowered, mirroring the lack of 
any information on radio objects at faint fluxes. We have 
chosen three models which span the range of reasonable 
N{z); two representing the extremes, and one intermediate. 

The first model we adopted is DP's 7 (called hereafter N3), 
indicated by Peacock (private communication) as the best 
model, which predicts N{z) for pure luminosity evolution. 
This is illustrated in figures |l| and ^ by the short-dashed 
line, and is characterised by a low-redshift component which 
becomes more dominant at fainter fluxes, and by a broad 
and shallow maximum for z > 0.5. The low-z spike could 
be interpreted as modelling the starburst-galaxy population 
which appears at these faint fluxes. 

The second model considered in our analysis, DP's 3, (called 
hereafter N2) is represented in figures |l| and |^ by the long- 
dashed line; it has a narrow and prominent bump around 
z ^ 1 and no low-redshift spike. 

As an intermediate model A''l we chose the average of DP's 
1-4, 6-7 for it shows both the low-redshift component and 
the 2: ~ 1 bump, but neither feature is too dominant. This 
is plotted as the solid line in Figures 1 and 2. 



3 PREDICTING THE ANGULAR 
CORRELATION FUNCTION 

The standard way of relating the angular two-point correla- 
tion function Wradioifi) to the spatial two-point correlation 
function £,radio{r, z) is by means of the relativistic Limber 
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Figure 2. Redshift distribution N(z) for the radio source pop- 
ulation at 1.4 GHz at 10 mjy. The dotted and dashed curves 
represent the 6 models (1-4, 6-7) taken from Dunlop & Peacock 
(1990); the solid curve is the average. 



Figure 3. The linear-theory density growth D{z) scaled to unity 
at the present time for Qq = 1, Hq = 0.3 and Qo = 0.1. The solid 
lines represent open models (A = 0) while the dashed lines are 
for flat models (A + = 1). 



equation (Peebles, 1980): 

Wradio(P) = ^ — — 72 , (J-j 

[f^ F-^{x)x^^{x)dx\ 

where x is the comoving coordinate, F{x) gives the correc- 
tion for curvature, and the selection function <^{x) is deter- 
mined by the N{z), 



N = I ^{x)F'^(x)x'^dx = 



N{z)dz, 



(2) 



in which A/" is the mean surface density on a surface of solid 
angle fla and N{z) is the number of objects in the given 
survey within the shell {z, z + dz). 

The physical separation between two sources separated by 
an angle 9 is given (in the small angle approximation) by: 



1 



2 \ 1/2 

(l + z) V +^ 



(3) 



The comoving coordinate x and the correction factor F{x) 
are different for different geometries; for a Universe with 
generic density parameter Qq and cosmological constant A = 
(see e.g. Treyer & Lahav, 1996): 



2c 
IF 

and 

F{x)^ 



floz - (Oo - 2)(1 - yr+T^i) 



ng(l + z) 



1 + ("0 



1/2 



(4) 



(5) 



This leads, with equation (□), to the following expression for 



the angular correlation function: 

oH0^2 N\z)/P{^l»,z)dz f^ ^rai^o{r,z)du 
Wradio(t^) — ^ "0 TT— ' 

{i::N{z)dz\ 

where P(Q.(i,z) is given by 

4(!^o - 1)[(1 + ^i^zfl'^ - 1] + n,^(l - z) + 2^02 



P(no,2) 
and u is defined by (B 



(l + z)2(l + J7o2)l/2 



In the case of a cosmological constant A 7^ with f2o + A = 1 
(flat space), we have F{x) = 1 and: 



x-—Qr^'^ 
no 



dz 



[(1-f 2)3+00^-1 



1/2 ' 



(8) 



(see Peebles, 1984; Treyer & Lahav, 1996) so that the ex- 
pression for wiO) assumes the form 

N'^{z)Q{^0,z)dz irad^o{r,z)du 

n...dU0)-2-n, [J-Niz)dzY 

with 

Q{Qo,z)^[{l + zf + no^-^'^^ (10) 



To take into account the evolution of clustering with epoch, 
we express the correlation function (.radioir, z) as: 

^radio{r,z) = ^ma3s{r)D'^{z)b^{z), (11) 

We have chosen the present-day ^massir) to have a shape 
which matches the observed ^gai from local optical galaxy 
surveys (eg. APM, LCRS, etc.) and a normalization which 
matches the 4-year COBE data (Bunn & White 1997). On 
large scales we use the linear-theory prediction from the 
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Figure 4. Integrated correlation function f(r, as a function of 
the angular scale 9 in the case Cq = 1 at redshift z = 0.1, 1, 2, 3, 4. 
Solid lines represent the predictions for growth of clustering under 
linear theory on the assumption of a i,mass (f") as described in the 
text, while the dashed lines take into account the effects of non- 
linear evolution. The two vertical lines bracket the range 0.3° < 
< 3° provided by the observations. 



power spectrum P{k) of Bond & Efstathiou (1984). On small 
scales the linear prediction under-estimates the true am- 
plitude, and so we extrapolated with a power-law of slope 
— 1.7, as observed from the APM correlation function (Mad- 
dox et al., 1990). For T — 0.2 the extrapolation is used for 



r < r, 



gal 



5.4Mpc. For other values of F we have re-scaled 



with a ■ 



^ so that 5i(r) oc ^2(fFl/F2). In each case we use the 4 
year COBE data to fix the normalization in terms of mass 
fluctuations a^"'' . 

D{z) is the linear-theory density growth-rate whose generic 
expression is given by Peebles (1980) and Lahav et al. (1991); 
D{z) = {1 + z)~^ for an Einstein-de Sitter Universe 

1 = (1 + z)~^ and f{z) ~ R°-^ , where 

R{z, fio. A) = !^o(l + zf[^o(l + zf - 

(no + A-l)(l + z)'+A]-' (13) 

Figure ^ shows the behaviour of D{z) for a range of Q.q 
and A. Our choice to consider only the case for growth 
of clustering under linear theory is extensively discussed 
in Magliocchetti et al., 1998; given the high median red- 
shift of objects in the FIRST survey ~ 1), the data 
obtained in the angular range 0.3° < 6* < 3° provided 
by our analysis of the catalogue correspond to a range of 
spatial scales (r > 10/i~^Mpc) where the linear regime 
still holds. As a check we calculated the integrated corre- 
lation function ^(r, z) = ^(y, z) dy as obtained for 
£,{r,z) = ^{r)D'^{z) (growth of clustering under linear the- 
ory) in the case of (,mass (r) as described earlier in this section 
and b{z) — 1. This integral was performed at different red- 



Figure 5. Trend of the quantity D{z) b{z) as a function of red- 
shift for three different models of biasing as described in the text. 
Qq is set to be 0.3, while b^^^ = 1.46. The solid lines represent 
open models while the dashed lines are for flat models. 



shifts and for f2 = 1 and the resulting integrated correlation 
functions have been compared with the expression 



X + 0.358a:^ + 0.0236a:^ 
1 + 0.0134x3 + 0.00202x9/2 



(14) 



from Hamilton et al., 1991 to take into account the effects 
of non-linear evolution. The results are shown in figure ^ 
where we plotted ^(r, z) as a function of the angular scale 
6 = r/xa with x given by equation ^ in the case flo = 1, for 
z = 0.1, 1, 2, 3, 4; solid lines represent the evolution of the in- 
tegrated correlation function under linear theory, while the 
dashed lines are obtained from equation ( |l^ . It is clear from 
the plot that differences between the two models occur only 
at angular scales significantly smaller than the range of our 
observations (0.3° < ^ < 3°) which are bracketed by the 
two dashed vertical lines. This gives some confidence that 
our model assumptions are reasonable. 

The biasing factor b{z) in equation ( |ll[ ) allows for evolution 
in the way radio sources trace the mass distribution. In prin- 
ciple b can be a function of both scale and redshift, but at 
z — 0, significant deviations from b(r) = const are seen only 
on small scales {r < 10/i"^Mpc) (Bagla, 1997). Thus in the 
following analysis we ignore any variation of bias with scale. 
We take into consideration three possible models for the 
evolution of biasing as a function of the redshift (for an ex- 
tensive study see Matarrese et al., 1997 and Moscardini et 
al., 1997). The first of these, called hereafter BO, is: 



b{z) = bo 



1.7/2 



(15) 



representing constant bias at all epochs. The two other mod- 
els are derived from the expression 



b{z) = fe-i + 



b-i){l + z) 



(16) 
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with 



(S) 



1.7/2 



Setting 6_i = /3 = 1 gives model Bl, while setting fe_i = 
0.41, l3 ~ 1.8 gives model B2. 

The Bl parameters have been chosen in agreement with 
the so-called galaxy conserving model (Fry, 1996; Tyson, 
1988) in which galaxies form at some particular redshift 
and then evolve without losing their identity, by following 
the continuity equation. The latter set (B2) refers to re- 
sults from N-body simulations (see e.g. Bagla, 1997) and 
describes the so-called merging model in which faint galaxies 
are subunits that merge to make up more luminous galaxies 
(Broadhurst et al., 1992; Clements & Couch, 1996; Baugh 
et al., 1996). ro is the clustering length of radio sources at 
z = (~ lOft-^Mpc in our case, as measured by Peacock 
& Nicholson, 1991; see also Magliocchetti et al., 1998) while 
Tq"' = 5.4/i~^Mpc is the corresponding value for optically- 
selected galaxies as measured in the APM survey, so that the 
factor ( j^) ^ allows for the different bias level of 'optical' 
galaxies compared to radio objects as seen at zero redshift. 
The mass r.m.s. fluctuation amplitude inside a sphere of 
8/i-'Mpc is given by a^""" , and the quantity b""' = 1/(7^""' 
is the bias of the distribution of optical galaxies relative to 
the distribution of mass. This value depends on both the ge- 
ometry and the normalisation of the power spectrum P(k). 
Note that the amplitude of £,mass{r) in equation (11) de- 
pends on the evolution of the biasing factor with redshift 
given in equations (15-16) through the relation: 



iradio{r,Q) = b^{fi)^,nass{r). 



(17) 



From equation (17) we can therefore determine the values 
for the mass correlation length r^^"-"" that turns out to be: 

1/7 

(18) 



mass 

ro 



where is a factor depending on the slope of the correlation 
function at small scales (see Peebles, 1980). 
As an illustration figure ^ shows the quantity D{z) b{z) as 
a function of redshift for the three biasing models described 
by equations (15-16) for Qo ~ 0.3, flat (dashed lines)/open 
(solid lines) Universe and 6q''' = 1.46. The evolution of 
clustering with redshift/epoch is completely determined by 
the combination of bias b{z) and density growth rate D{z) 
(equation ^l|), and so the different forms of b{z), give models 
where a) clustering decreasing with look-back time (BO) , b) 
clustering roughly constant (51) and c) clustering rapidly 
increasing with look-back time (-B2). 

To constrain an already wide parameter space we flx the 
value for the reduced Ifubble constant as ho = 0.65 (see 
e.g. Kundic et al., 1997). Moreover, it is almost impossible 
to discriminate between different values of fio through the 
evolution of clustering, because low density models in which 
clustering evolves slowly have a geometry which provides 
more time for evolution. This degeneracy in flo is often re- 
ferred to as cosmic conspiracy. Hence we restrict our analysis 
to r^o ~ 0.4, allowing also for flat geometries {Qo -|- A = 1). 
We adopt the normalisation of the r.m.s. mass fluctuations, 
0-™"°", determined from the 4-year COBE data by Bunn & 
White (1997). So choosing a well-deflned cosmology (i.e. a 



fixed value of fio and A) fixes the value of ag and therefore 
the amplitude of £_mass{r), given by r™'"'^ from equation 
(18) Then the parameter b'g'* in the expression for b{z) is 
also determined by equation (15) and (16). We then have: 



Qo = 0.4 A = 
flo = 0.4 A = 0.6 



mass 



mass 



0.64; rj^ 
1.07; r^ 



3.31/i"^Mpc 
6.23/i"^Mpc 



As a comparison for a standard CDM model we would have 
^n^ass ^ -^ 22; rj""'" = 6.82/i-iMpc. Thus it turns out that 
in ACDM and standard CDM (SCDM) models the distribu- 
tion of optical galaxies is "antibiased" relative to the distri- 
bution of mass so that the amplitude of the corresponding 
^massj-^-j -g jjjgjjgj. than that found for galaxies; the opposite 
case pertains for open CDM (OCDM) models. 
Concerning the constant bias case BO we have instead de- 
cided to fix (ad hoc) aj""'" = 1, so that r^^""" = 5.4 h'^Mpc. 

Note that in the case of F = 0.5 both r^"-"" and cr^""" wiU 
differ from the values expressed for F = 0.2. In more detail 
we have: 

fio = 0.4 A = ^ o-^""" = 1.51; = 8.95/i''Mpc 

fio = 0.4 A = 0.6 (7^"" = 2.52; r^""" = 15.8/i~'Mpc, 

i.e. ff F = 0.5, in both the ACDM and OCDM models the 
distribution of galaxies is "antibiased" relative to the distri- 
bution of mass. 



4 PREDICTIONS FOR W{e) 

In this section we will focus our attention on the predicted 
angular correlation function w{0) obtained for two different 
CDM models with respectively, F = 0.5 (standard CDM 
model) and F = 0.2 (modified CDM model), values that 
bracket the range of models providing the best fit to the 
present observations of the power spectrum P{k) (Peacock 
& Dodds, 1994). Note that these models treat F as a free 
parameter independent of hoQo, and so do not represent 
consistent CDM models. 

We will focus on the three different redshift distributions of 
radio sources (Nl, N2, N3) described in section 2, and on 
the three functional forms for the evolution of the biasing 
with redshift (BO, Bl, B2) introduced in section 3. In order 
to test the self-consistency of the predictions of each model, 
we consider two different flux-cuts: 3mJy, the lowest flux 
limit where the incompleteness of the survey is negligible 
(see Magliocchetti et al., 1998) - and lOmJy, the brightest 
flux limit for which the measurements of Wobs{9) are not 
dominated by the errors). In what follows we always start 
by considering models derived for F > lOmJy, because the 
differences between the different N(z)^s are less significant 
than for fainter fluxes (see figures ^ and |^. In particular 
the relative weight of the low-redshift component and of the 
2 ~ 1 bump are much less variable at higher fiux limits. 

We will first show the results for a CDM power spectrum 
with F = 0.2 and then compare them to those with F = 0.5. 
The theoretical estimates of w{9) obtained at lOmJy are 
to be compared with those for 3mJy, while keeping all the 
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Figure 6. Theoretical prediction of the angular correlation function for T = 0.2, N{z) s ATI, N2, N3 at lOmJy and 3mjy, and different 
bias models as described in the text. The solid lines are for open models while the dashed lines indicate spatially-flat models. 



other parameters fixed, in order to test tlie consistency of 
the predictions of different models for N{z). 

In Figure 6 wo plot the theoretical estimates of w{6) (in 
the case F = 0.2) at lOmJy (upper panel) and 3mJy (lower 
panel), for the three models of N{z) introduced in section 2 
and the three models for the evolution of bias with redshift 
of section 3. The correlation length ro in equations 15 and 
16 has been fixed at 10h~^ Mpc in agreement with both 
our previous estimate (Magliocchetti et al., 1998) and the 
results found by Peacock and Nicholson (1991). The solid 
lines are for open (f2 = 0.4) models while the dashed lines 
are for flat (A ^ 0) geometries. 

Independent of the flux cut, the discrimination amongst dif- 
ferent models comes both from the eventual fall of w(6) at 
large angular scales and from its overall amplitude. The 
overall amplitude of the angular correlation function is 
mainly related to the level of bias at high z: models char- 



acterised by a bias factor increasing with redshift (B1/B2) 
show a stronger clustering signal, especially at smaller angu- 
lar scales, relative to models with constant bias (BO). Their 
amplitude increases as the dependence of b{z) on the factor 
(1 -|- z) in equation (16) becomes steeper, as can be seen by 
comparing the B2 with the Bl models. The difference in am- 
plitude between the BO and Bl models is not pronounced; 
this effect is primarily due to the lower normalisation of 
the Bl models deriving from having taken a^"'"" 7^ 1 as 
opposed to the assumption of = 1 in the BO case. 

Furthermore, as the main contribution to the clustering sig- 
nal at large angular scales is given by low-redshift objects 
(see Cress & Kamionkowsky, 1998), the factor (1 + z) is not 
strong enough to drive the evolution of the Bl models far 
from that described by BO. 

The fall at large angular scales depends instead on the neg- 
ative feature in ^, a genuine feature of CDM-like models, 
as well as on the amplitude of the low-z component in the 
N{z). Note that in equations 6 and 9 the integral is weighted 
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Figure 7. Theoretical prediction of the angular correlation function for T = 0.5, N{z) = Nl, N2, N3 at lOmJy and 3mjy, and different 
bias models as shown in the figure. The solid lines are for open models while the dashed lines indicate spatially flat-models. 



by N{z)'^, so the spike makes a significant contribution to 
w{6) even though the volume is small. The low-z amplitude 
describes the relative contribution of the local population of 
radio objects, that, as already niciitioiicd, is the component 
which provides the dominant contribution to the cluster- 
ing signal at large angular scales (Cress & Kamionkowsky, 
1998). In fact from figure 6 one can see how all the correla- 
tion functions obtained for N2 (no low-redshift component) 
drop at 6^ ~ 2°, while those for N3 and Nl models, charac- 
terised by a larger low-redshift component, flatten out at this 
angular scale, with a positive tail at larger scales that be- 
comes more predominant as the number of nearby objects is 
increased. This effect is stronger at 3mJy than at lOmJy be- 
cause the amplitude of the local population of star-bursting 
galaxies increases as the flux threshold is lowered. 

The predictions for w{9) for the case of F = 0.5 are illus- 
trated in figure 7; the models show basically the same fea^ 
tures as those obtained in the analysis for F = 0.2, indicating 



very little dependence on the value of the shape parameter 
F. 



5 RESULTS FROM THE COUNTS IN CELLS 
ANALYSIS AT DIFFERENT FLUX LIMITS 

This section considers the counts-in-cells analysis of the 
FIRST catalogue carried out at different flux limits; we then 
compare the predictions with the corresponding measure- 
ments in section 6. We measured the clustering of radio ob- 
jects using the so-called galaxy distribution function that 
gives the probability of finding A'' objects in a cell of par- 
ticular size and shape. By defining the fe-th moment of the 
counts as 

Mfc = ((iV-iV)'=) (19) 
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Figure 8. The normalised variance cr^ vs the cell size © for the FIRST survey for objects brighter than 5mjy (left panel), 7mjy (middle 
panel) and 10m Jy (right panel). Errors are estimated from the variance in four random subsets. 
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Figure 9. Amplitude of the angular correlation function w^bsiS) 
as measured from the FIRST survey as a function of the flux 
limit. 



where N = nQ. is the mean count in the solid angle Q,, we 
have that the second moment /i2 of the galaxy distribution 
function is related to the two-point correlation function w\2 
through the expression 




where is the shot noise resulting from the discrete nature 
of the sources (Poisson noise), and 



*2 



^2 



w-12 dQi dQ2 



(21) 



is the normalised variance in terms of the two-point correla- 
tion function integrated over a cell of area and particular 
shape. By assuming a power-law form 



w{e) = A 6^' 



(22) 



and by considering square cells of size 57 = O x O square 
degrees, we can evaluate the integral in equation (j^) (see 
Totsuji & Kihara, 1969), obtaining 



^2 



TV 



Ai 



''dQ.^dQ.2 



ACQ" 



(23) 



^12 = N + (Ny^2 



(20) 



where 6*7(7) is a coefficient which can be evaluated numer- 
ically by Monte Carlo methods (Lahav & Saslaw, 1992). It 
is therefore possible to use the — Q relation to evaluate 
the two parameters {A, 7) which describe the correlation 
function (22). 

The following analysis has been carried out as in Maglioc- 
chetti et al. (1998) for different versions of the original 
FIRST catalogue obtained by combining source components 
into single source, following well defined criteria. The results 
are shown in figure ^ where we plot as a function of Q 
for difi'erent flux limits. The slopes and the intercepts of the 
plots are estimated by a least-squares procedure minimising 
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the quantity 

X-(a,5) = f (24) 

i=l ^ ' ^ 

with b = (5 — 7), a = \og{AC.y). The errors Ai are obtained 
using the 'partition bootstrap method' in which the nor- 
malised variance is calculated for four subdivisions of the 
survey region and the standard deviation of these measure- 
ments at each angle is used as a measure of the error. From 
this analysis we find, for F > 3mJy, F > 5mJy, F > 7mJy 
and F > lOmJy respectively: 



7 


= 2.50 ±0.1 ; 


A = 


(1.06 ±0.1) • 


10~^ 


(25) 


7 


= 2.50 ±0.1 ; 


A = 


(1.00 ±0.1) ■ 


10~^ 


(26) 


7 


= 2.50 ±0.09 


; A-- 


= (0.97 ±0.1) 


• lo"^ 


(27) 


7 


= 2.40 ±0.2 ; 


A = 


(0.89 ±0.1) ■ 


lo"^ 


(28) 



The 3mJy measurement is taken from Magliocchetti et al. 
(1998). Above 10m Jy the measurements are dominated by 
Poisson errors. 

From equations (Q, (|^ and ( |28[ ) we have that the slope 7 
of the correlation function (see eq. is independent of the 
flux limit in the range 3mJy< F < lOmJy (see Magliocchetti 
et al., 1998 for the results at SmJy). Furthermore we find 
that the amplitudes A at each flux-cut are vary smoothly 
with the flux limit, as shown in figure ^ We can therefore 
be confident that the clustering signal found in our analysis 
is mainly given by populations of radio objects which do 
not change their clustering properties as the flux level cut 
is increased. This result will be of particular help in the 
next sections when we compare the theoretical predictions 
for the correlation function with measurements at different 
flux limits. 



6 COMPARISONS WITH THE DATA 

One of the main reasons for preferring the variance statistic 
(i.e. the counts in cells analysis) to the the direct evaluation 
of the correlation function S^radio is that is a cumulative 
quantity which provides a complete description of large-scale 
structure in a more robust way. Therefore, in order to match 
the models obtained for the angular correlation function in 
section 4 with the results inferred from the counts in cells 
analysis, one has to integrate w{9) to get the predictions 
for the variance of the distribution of sources. This can be 
done by means of equation (23) which, for a generic w{9), is 
written as 

cr^ = y" w{e) SQi 5^2, (29) 

and in the case of square cells, can be expressed as the two 
dimensional integral 

cr^(e) = eM dx w{e)dy, (30) 

Jo Jo 

where 0^ is the area of the cells and 9 = ^yx^ + y^. The 
quantity given by equation (29) has then been evaluated 



for all the w(9) models illustrated in section 4 and the cor- 
responding (T^'s have been compared with the estimates of 
(Toftj obtained from the data both at 3mJy and lOmJy as 
explained in section 5 and in Magliocchetti et al. (1998). 

Figures 10 and 11 show the results for F = 0.2 and F = 0.5 
respectively; the symbols are as in section 4 and once again 
the solid lines are for open models while the dashed lines are 
for spatially-flat models. Independent of the functional form 
adopted for the redshift distribution A'^(2) and of the flux 
cut, models with bias strongly evolving with redshift can be 
ruled by the data with a high confidence level (> 3cr) as they 
grossly overestimate the amplitude of u^{0) at all scales. 
Although the predictions for constant bias (BO) provide a 
somewhat better fit to the data (particularly for F = 0.2), 
the BO (constant bias) and the Bl (bias linearly evolving 
with redshift) models are almost degenerate with reference 
to the measurements. As already mentioned in section 4, 
this effect is due to the low normalisation of the mass cor- 
relation function S^massir) in the Bl case, resulting from the 
assumption that optical galaxies are biased relative to the 
distribution of mass, and from the fact that for Bl (contrary 
to the B2 models) the dependence of b{z) on the redshift is 
not strong enough to drive the evolution of ^radioir,z) to 
higher amplitudes, on these angular scales. 
If one assumes that radio sources are indeed strongly biased 
with respect to "normal" galaxies by the factor (10/5.4)^'^ 
appearing in equations (15) and (16), then the conclusion of 
the BO model as favoured for the interpretation of the data 
is in agreement with the results found by Matarrese at al. 
(1997) and Frieman & Gaztanaga (1994) from the angular 
correlation function in optical/infrared surveys. 

This kind of analysis apparently is not sensitive to differ- 
ences in N{z)'s. As ct^ is an integral quantity, its behaviour 
has very little dependence on the form of the tail of ■w{9) 
at large angular scales where the clustering signal is weaker 
by at least two orders of magnitude than that at S ~ 0.1°. 
Thus, although the models for w{9) for different redshift 
distributions show remarkable differences, the combined ef- 
fects of the integral (28) and the presence of observational 
uncertainties in the estimates of o-^,j^(0) wash out these dif- 
ferences. 



7 THE CLUSTERING EVOLUTION FOR 
DIFFERENT RADIO POPULATIONS 

A major problem in the analysis and interpretation of the 
clustering of radio sources at faint flux limits is the mix- 
ing of different populations. The brighter objects in a deep 
radio survey are identified with powerful radio-loud QSOs 
and AGN's; the fainter objects (< 3mJy) are dominated by 
a population of faint low-z sources which are local starburst- 
ing galaxies (see Wall, 1994). As radio-loud QSOs tend to 
be hosted by giant ellipticals within rich clusters (see e.g. 
Cristiani, 1998), these objects show a much stronger cor- 
relation than normal field spirals (see also Loveday et al., 
1995) containing the population of star-bursting galaxies. 
It is therefore of crucial importance to take into account 
the characteristics of these two different populations to get 
information on their clustering properties and on the evo- 
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Figure 10. Theoretical predictions for the normalised variance in the case of T = 0.2, N{z) = ATI, N2, N3 at lOmJy and 3mjy, and 
different bias models. The solid lines are for open models while the dashed lines are for spatially-flat models. The data are from the 
FIRST survey 



lution of bias with redshift. Theoretical works suggest that 
bias depends on the mass of the dark matter halos in which 

galaxies were born (Mo & White, 1996) and therefore its 
dependence can also be related to the mass of the galaxy it- 
self in the hypothesis that the most "visible" galaxies reside 
within the most massive dark matter halos. 

To perform our analysis we use a simple approach assum- 
ing at low redshifts the clustering signal is dominated by the 
starbursting objects while at higher z's both the populations 
contribute. This approximation is justified by the observa- 
tional evidence of a drastic drop in the number of quasars 
in the nearby universe. We therefore write the correlation 
function as follows: 

itot{r, z) = {A'U + 2AB(5aSs)''" + B'iB ] (31) 

where E,A(r, z) is the correlation function for the starbursting 
galaxies, ^^(f, z) is for radio- loud quasars, and A and B arc 
the relative space densities the two populations normalized 
so that A + B = 1. We assume that B = for 2 < 0.1, and 
consider three possibilities for z > 0.1: starburst dominated. 



A = 0.8, B = 0.2; an equal mix A = 0.5, B = 0.5; and AGN 
dominated, A = 0.2, B = 0.8. For the starburst correlation 
function wc use ^A{r,z) — ^ D{z)^b' (z)'^ , where ^(r) and 
D{z) are as expressed in section 3. For the AGN correlation 
function we use ^B(r, 2) = D{z)^b{z)^ where ^'(r) is the 
effective linear bias mass correlation function for QSOs: we 
assume on large scales that it is the linear prediction for a 
CDM model with generic F, while at small scales ^ has been 
extrapolated as a power-law of slope —2.1 as found from 
the analysis of the clustering of bright-early type galaxies 
in the Stromlo-APM redshift survey (Loveday et al., 1995). 
Note that the value of 2.1 for the slope is different from 
that used in section 3, where we considered the radio sam- 
ple to be homogeneous. b{z) and b'{z) are the bias associ- 
ated with bright radio sources (QSOs -I- AGN) and spiral 
galaxies (starbursting population) relative to the distribu- 
tion of mass; their expressions are given by equations (15) 
and (16) assuming ro ~ lOft^^Mpc (powerful radio sources) 
and ro ~ 5.4/i~^Mpc (optical-galaxies). This yields for ex- 
ample, in the case of constant bias (BO) models, that the 
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Figure 11. Theoretical predictions for the normahsed variance cr^ in the case of F = 0.5, N{z) = Nl, N2, N3 at lOmJy and 3mjy, and 
different bias models. The solid lines are for open models while the dashed lines are for spatially-flat models. The data are from the 
FIRST survey. 



starbursting population is indeed unbiased (6(«) = 1), while 

the AGN population shows a higher bias relative to the dis- 
tribution of the mass; eX z = 0, b{0) = {-^Y'^^^ -^Jaww- The 

expression for ^tot{r,z) has then been projected by means 

of the Limber equation (1) in order to got the relative func- 
tional form for the angular correlation function w{d). 

As Cress and Kamionkowsky (1998) have indicated, at scales 
S > 1° the angular correlation function is dominated by 
nearby objects (within z < 0.1), while at scales 8 ~ 0.1° 
there is an equal contribution from local {z < 0.1) and dis- 
tant {z > 0.1) sources. Thus, to study the evolution of the 
clustering of bright distant objects we have to compare our 
models of w(0) with the observed correlation function push- 
ing the measurements down to angles 9 ~ 0.02°. Note that, 
as can be seen in figure 4, the assumption of growth of clus- 
tering under linear theory is still a good approximation at 



these angular scales. The counts-in-cells statistics for the 
FIRST survey are dominated by shot (Poisson) noise for 
9 < 0.3° (see Magliocchetti et al., 1998) and so to extend 
the range of measurement down to smaller scales we have 
measured angular correlation function Wobs{9) directly from 
the catalogue. 

Wc recall here that the angular correlation function wi2 — 
w(9) gives the excess of probability, with respect to a Pois- 
son random distribution, of finding two sources in the solid 
angles 5Qi 5^2 separated by an angle 9, and it is defined as 

6P = n'^dili6n2 [1 + w{9)] (32) 

where n is the mean number density of objects in the cata- 
logue under consideration. We measured Wobs from the esti- 
mator (Hamilton, 1993) 

Wobs{9)=4^^^ -1 (33) 
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Figure 12. Theoretical prediction for the angular correlation function in the case of F = 0.2, N{z) = Nl, N2, N3 at 3mjy, and different 
bias models as shown in the figure. The solid lines are for models with A = B = 0.5, the long-dashed lines are for A = 0.8, B = 0.2, 
while the short-dashed lines stand for A = 0.2 and B = 0.8 (see text for details). The data are obtained from measurements of Wabs{S) 
using equation pA. The error bars show the Poisson estimates for the catalogue under consideration. 



where DD is the number of distinct data-data pairs, DR 
is the number of random-data pairs and RR is the num- 
ber of random-random pairs as a function of angular sep- 
aration and where the random data have been obtained 
by generating a random catalogue of positions and selected 
sources above the sensitivity limit from the coverage map 
(see Magliocchetti et al., 1998 for further information about 
the catalogue). 

Given the very weak dependence of the angular correlation 
function on geometry, the shape parameter of the power 
spectrum F and the flux cut (at least for F < lOmJy), we 
restrict our analysis ton = l,r = 0.2 and F > 3m Jy, while 
still allowing for different functional forms of the redshift 
distribution N{z). All the models assume that the distri- 
bution of starburst objects is unbiased with respect to the 
mass distribution, justified by the results obtained for the 
more local population in section 6. The corresponding mod- 
els are shown in figure 12 where we also plot Wobs taken 
from the FIRST survey at 3mJy; the error bars show Pois- 
son estimates for the catalogue under consideration. The 
solid lines are for models in which the distribution of the 
AGN population is unbiased (BO), has a bias linearly evolv- 
ing with redshift (Bl) and has a bias strongly evolving with 
z (B2). The predictions have been calculated in the case of 
the two populations equally contributing to the clustering 
signal {A — B = 0.5). As figure 12 shows, while at large 
angular scales the BO and Bl models are roughly degener- 
ate as in section 6, pushing the analysis down to smaller 
scales removes this degeneracy. In fact the BO and B2 mod- 



els are ruled out as possible descriptions of the data as they 
show either too small (BO) or too large (B2) amplitude, es- 
pecially at intermediate scales. The best fit is provided by 
the Bl curves i.e. by models in which the distribution of 
the starbursting population is taken to be roughly unbiased 
while that one for the AGN population shows a bias linearly 
evolving with redshift. 

The effect of varying the relative weight of the two pop- 
ulations in the Bl model is shown in figure 12 for the two 
cases: ~ 20% starburst and ~ 80% AGN (short dashed lines) 
and ~ 80% starburst and ~ 20% AGN (long dashed fines). 
The agreement with the data is best for the AGN dominated 
model. The notable differences for different N{zYs found in 
section 4 reside primarily in the part of the angular correla- 
tion function at large scales where the signal is dominated 
by the noise. 

The fact that we need a biased distribution of AGN in or- 
der to match the data is in agreement with the hierarchical 
models of galaxy formation predicting the less massive ob- 
jects to be less strongly clustered than those formed in the 
higher and biased peaks of the density field (Giavalisco et 
al., 1998), with the most powerful sources being hosted by 
the most massive dark matter halos, and the biasing fac- 
tor depending (at least) both on the epoch and the type of 
objects under consideration (i.e. on the halo mass). 
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8 CONCLUSIONS 

To put constraints on the redshift distribution N{z) of radio 
objects at faint fluxes and on their clustering and biasing 
properties, we have calculated predictions for the angular 
correlation function w{9) using three models for N{z) and 
three models for the clustering evolution. The three models 
for the redshift distribution have been chosen to span the 
range of reasonable possibilities. The predictions for w{9) 
have been calculated for a ^mass obtained from two CDM 
models (F — 0.2 and F = 0.5), three different assumptions 
for the evolution of the bias and in the case Qq ~ 0.4, both 
flat and open space. We chose a fixed value for the density 
parameter because of the degeneracy amongst predictions 
obtained for different Qq's. The bias models correspond to 
three sensible galaxy evolution models (galaxies tracing the 
mass at all epochs - constant bias; b{z) oc (f + z) - galaxy 
conserving model; b{z) oc + * - galaxy merging model). 

The discrimination amongst different models comes both 
from the fall of the calculated w{6) at large angular scales 
and from its overall amplitude. The fall depends on the neg- 
ative feature in ^ at large scales, which is a generic feature 
of CDM-like models, and also on the low-z component in 
the N(z) whose amplitude depends on the presence and rel- 
ative contribution of any local population of radio objects. 
The overall amplitude of the angular correlation function is 
instead mainly related to the level of bias at high redshift. 



sources, with clustering properties very similar to normal 
spiral galaxies and a distribution that is roughly unbiased 
relative to the distribution of mass, and b) a more numer- 
ous population of radio- loud AGN's much more strongly 
clustered and whose distribution is biased relative to the 
distribution of the mass, with a bias factor which evolves 
linearly with redshift. This result is in good agreement with 
hierarchical CDM theories of galaxy formation that predict 
the lower luminosity /lower mass objects being less clustered 
than the high- luminosity sources, with the brighter sources 
being associated with higher (and biased) peaks in the dark 
matter distribution. 

There are crucial observations needed to further the analy- 
sis: in particular observational measurements of N(z), from 
identifications and redshift measurements for a complete 
sample, would be of great importance in constraining struc- 
ture formation models as well as understanding the distri- 
bution and evolution with time of the different populations 
of radio objects at mjy levels. 
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For a comparison with the data inferred from the counts-in- 
cells analysis of the FIRST survey we integrated the models 
for w{8) to calculate the corresponding predictions for the 
variance of the distribution ; we then matched the results 
with the measurements for different flux limits presented in 
section 1 and in Magliocchetti et al. (1998). Although mod- 
els for the angular correlation function obtained for different 
N{z) show quite remarkable differences, in the framework of 
this analysis it is not possible to discriminate amongst them, 
given the errors in the measurements and the fact that the 
integral used to calculate washes out the differences. 
Concerning the bias evolution, all the models of bias strongly 
dependent on redshift can be ruled out with high confidence 
level as they grossly overestimate the amplitude of at 
all scales. There is a degeneracy between models with con- 
stant bias and those for bias depending linearly on redshift, 
although the predictions for constant bias are mildly pre- 
ferred by the data, especially for F = 0.2. The analysis is 
quite insensitive to the geometry of the Universe. 

As a further step, we introduced a model for the correla- 
tion function able to account for a mixing of populations 
with completely different properties, the first population of 
faint starbursting/spiral galaxies, the second one of AGN- 
powered objects normally hosted by bright ellipticals. We 
assumed that starbursting population traces the distribu- 
tion of mass at all times, and we allowed for biasing evolu- 
tion in the case of radio-loud AGN. Appropriate measure- 
ments for Wobsid) were made directly from the catalogue. 



REFERENCES 



It turn ) out that an evolution ot bias with redshift is re- 
quired tor aGN to match the data, ihe model that best 
fits the observed angular correlation function is given by the 
sum of a) a less numerous population of faint starbursting 



Bagla, J.S., 1997; |astro-ph/9711081 
Baugh CM., Cole S., Frenk C.S., 1996; MNRAS, 282, L27 
Bond J.R., Efstathiou G., 1984; ApJ, 285, L45 
Becker R.H., White R.L., Helfand D.J., 1995; ApJ, 450, 559 
Broadhurst T.J., Ellis R.S., Glazebrook K., 1992; Nat, 355, 55 
Bunn E.F., White M., 1997; ApJ, 480, 6 
Clements D.L., Couch W.J., 1996; MNRAS, 280, L43 
Cress CM., Helfand D.J., Becker R.H., Gregg M.D., White R.L., 

1996; ApJ, 473, 7 
Cress CM., Kamionkowsky M., 1998; istro-ph/9801284 
Cristiani S., 1998; 1998elss.con£E..41C 
Dunlop J.S., Peacock J. A., 1990; MNRAS, 247, 19 
Frieman J. A., Gaztaiiaga E., 1994; ApJ, 425, 392 
Fry J.N. 1996; ApJ, 461, L65 
Giavalisco M. et al. 1998b; in preparation 
Hamilton A.J.S., 1993; ApJ, 417, 19 

Hamilton A.J.S., Kumar P., Lu E., Matthews A. 1991; ApJ, 374, 
LI 

Kooiman L.K., Burns J.O., Klypin A.A., 1995; ApJ, 448, 500 
Kundic T., Turner E.L., CoUey W.N., Gott J.R.IH, Rhoads J.E., 

Wang Y., Bergeron L.E., Gloria K.A., Long D.C, Malhotra 

S., Wambsganss J., 1997; ApJ, 482, 75 
Lahav O., Liljc P.B., Primack J.R., Rees M.J., 1991; MNRAS, 

251, 128 

Lahav O., Saslaw W.C, 1992; ApJ, 396, 430 
Loan A. J., Wall J.V., Lahav O., 1997; MNRAS, 286, 994 
Loveday J., Maddox S.J., Egstathiou G., Peterson B.A, 
ApJ, 442, 457 

Maddox S.J., Efstathiou C, Sutherland W.J., Loveday J. 

MNRAS, 242, 43P 
Mag-Hocchetti M .. Maddox S.J., Lahav O., Wall J.V., 1998; kstro 



1995 



1990; 



ph/9802269 



Matarrcse S., Coles P., Lucchin F., Moscardini L., 1997; MNRAS, 
286, 115 

Mo H.J., White S.D.M., 1996; MNRAS, 282, 347 



© 0000 RAS, MNRAS 000, 000-000 



14 M. MagliocchetU, S.J. Maddox, 0. Lahav, J.V. Wall 



Moscardini L.. Coles P., Lucchin F., Matarrese S., 1997; astro- 



ph/^712184 

Peacock J. A., Nicholson D., 1991; MNRAS, 253, 307 
Peacock J. A., Dodds S.J. 1994; MNRAS, 267, 1020 
Peacock J. A., 1997; MNRAS, 284, 885 

Peebles P.J.E., 1980; The Large-Scale Structure of the Universe, 

Princeton University Press 
Peebles P.J.E., 1984; ApJ, 284, 439 
Seldner M., Peebles P.J.E., 1981; MNRAS, 194, 251 
Shaver P.A., Pierre M., 1989; A&A, 220, 35 
Sicotte H, 1995; PhD thesis, Princeton University 
Steidel C.C., Giavalisco M., Pettini M., Dickinson M., Adelberger 

K.L., 1996; ApJ, 462, L17 
Steidel C.C., Adelberger K.L . Dickinson M.. G iavalisco M., 



Pet- 



atro-ph/9708125 



Giavalisco M.. Dickinson M., Pet- 



atro-ph/9805267 



tini M., Kellog M., 1997; 
Steidel C.C., Adelberger K.L 

tini M., Kellog M., 1998; 
Tegmark M., Peebles P.J.E., 1998; astro-ph/980407 
Totsuji H., Kihara T., 1969; PASJ, 21, 221 
Treyer M.A., Lahav O. 1996; MNRAS, 280, 469 
Tyson J. A., 1988; ^J, 91, 1 
Wall J.V., 1994; Austr. J. Phys., 47, 625 
Wall J.V., Jackson C.A., 1997; MNRAS, 290 L17 
Windhorst R.A., Miley G.K., Owen F.N., Kron R.G., 

1985; ApJ, 289, 494 



Koo R.C., 



© 0000 RAS, MNRAS 000, 000-000 



